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ABSTRACT 

High redshift galaxy clusters allow us to examine galaxy formation in extreme environments. 
Here we compile data for 15 2 > 1 galaxy clusters to test the predictions from a state-of-the- 
art semi-analytical model of galaxy formation. The model gives a good match to the slope and 
zero-point of the cluster red sequence. The model is able to match the cluster galaxy luminos¬ 
ity function at faint and bright magnitudes, but under-estimates the number of galaxies around 
the break in the cluster luminosity function. We find that simply assuming a weaker dust atten¬ 
uation improves the model predictions for the cluster galaxy luminosity function, but worsens 
the predictions for the red sequence at bright magnitudes. Examination of the properties of the 
bright cluster galaxies suggests that the default dust attenuation is large due to these galaxies 
having large reservoirs of cold gas as well as small radii. We find that matching the luminosity 
function and colours of high redshift cluster galaxies, whilst remaining consistent with local 
observations, poses a challenge for galaxy formation models. 

Key words: galaxies: abundances; galaxies: clusters: general; galaxies: evolution; galaxies: 
high-redshift; galaxies: luminosity function; methods: numerical 


1 INTRODUCTION 

Galaxy clusters are the most massive bound structures found in the 
Universe. Not only are clusters excellent proxies for massive dark 
matter halos (and therefore a useful cosmological probe), but they 
are also unique sites of galaxy evolution (e.g. Kravtsov & Borgani 
2012 ) 

Over the past decade a wealth of observational data has been 
gathered that points towards the redshift range z ~ 1.5 — 2.5 as be¬ 
ing a pivotal epoch in the evolutionary history of the galaxy popu¬ 
lation, with star formation, black hole accretion and galaxy mergers 
reaching their peak activity before being subsequently suppressed 
(e.g. Dickinson et al. 2003; Hopkins 2004). Similarly, observations 
indicate that at this epoch galaxy clusters and proto-clusters were 
in the process of being transformed from dynamical, merger-driven 
over-densities to the more relaxed systems that we see today. As 
such, there is mounting evidence that this redshift range marks the 
quenching of star formation in massive cluster galaxies and the 
build-up of the cluster red sequence (RS, Bower et al. 1992a,b; 
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Lidman et al. 2008; Brammer et al. 2011). However, developing a 
physical description of the rapid evolution of the galaxy population 
in clusters in this transformation phase represents a challenge for 
current models of galaxy formation. 

The next generation of cosmological galaxy surveys, such as 
the Dark Energy Survey (DES, The Dark Energy Survey Collabo¬ 
ration 2005) or the European Space Agency’s Euclid mission (Lau- 
reijs et al. 2011), are expected to observe many thousands of high 
redshift galaxy clusters. However, identifying galaxy clusters in the 
huge volumes probed by these surveys is a difficult task fraught 
with systematics, especially as these surveys will be predominantly 
photometric. The large uncertainties inherent in photometric red¬ 
shift estimation make the identification of cluster members based 
upon their spatial separation, as is the case for th q friends-of-friends 
(Huchra & Geller 1982) and Voroni-Delaney methods (Marinoni 
et al. 2002), much more challenging (e.g. Zandivarez et al. 2014). 

A more favourable approach for photometric datasets is to 
identify cluster members based upon whether they lie on the clus¬ 
ter RS in the colour-magnitude relation (CMR). Under the assump¬ 
tion that early-type galaxies dominate the cluster galaxy population 
and that this population follows a tight relation in colour-magnitude 
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space, then, when imaged in two photometric bands bracketing the 
4000A break, the cluster galaxies will be the brightest, reddest ob¬ 
jects (Stanford et al. 1998; Gladders & Yee 2000). The observa¬ 
tional efficiency of this approach has led to it being used exten¬ 
sively in cluster detection (Gladders & Yee 2005; Wilson et al. 
2006) and adopted in several group and cluster finding algorithms 
(e.g. Koester et al. 2007; Murphy et al. 2012; Rykoff et al. 2014). 
As discussed in Gladders & Yee, optical and infrared imaging of 
local and z > 1 galaxy clusters indicates the universal presence of 
a RS (e.g. Baldry et al. 2004; Bell et al. 2004; Brinchmann et al. 
2004), with many studies supporting a high formation redshift of 
the stellar population of Zf > 2 (e.g. Ellis et al. 1997; Smail et al. 
1998; Stanford et al. 1998; Ponman et al. 1999; Lopez-Cruz et al. 
2004; Gladders & Yee 2005; Miller et al. 2005; Voit 2005; Mei 
et al. 2006a,b; Koester et al. 2007; Gilbank et al. 2008; Lidman 
et al. 2008; Mei et al. 2009; Wilson et al. 2009; Gilbank et al. 2011; 
Lin et al. 2012; Mei et al. 2012). There is, however, some obser¬ 
vational evidence for ongoing star-formation in clusters at z > 1, 
(e.g. Hayashi et al. 2010; Brodwin et al. 2013; Fassbender et al. 
2014; Mei et al. 2015). 

Although using galaxy colours to identify galaxy clusters is 
preferable when dealing with the larger uncertainties inherent in 
photometric redshifts (e.g. Abdalla et al. 2011), this approach is 
sensitive to our understanding of the astrophysical processes gov¬ 
erning the evolution of galaxies in cluster environments and the 
build-up of the RS, as well as possible biases introduced by photo¬ 
metric colour selections. As a result, the RS method must be tested 
and calibrated. At low redshift this can be done with spectroscopic 
datasets, but at higher redshifts, where the spectroscopic data is 
sparse, one must turn to using synthetic ‘mock’ catalogues based 
upon the latest galaxy formation models (Baugh 2008). 

Due to their rarity, generating a population of galaxy clus¬ 
ters requires very large volume cosmological N-body simulations, 
which due to limitations in computational resources are typically 
dark matter only simulations. Several techniques are available for 
populating the halos from N-body simulations with galaxies. A 
common approach is to use empirical methods, such as the halo 
occupation distribution (HOD, Berlind & Weinberg 2002) or sub¬ 
halo abundance matching (Vale & Ostriker 2004). These methods 
have the benefit that they are tuned using observations to ensure 
that the luminosity function and colour distribution of the galaxies 
are correct by construction. At high redshifts, however, the lack of 
observational data prohibits the use of such methods, though some 
redshift dependent approaches have been proposed (e.g. Moster 
et al. 2013). 

Instead, semi-analytic galaxy formation models provide a 
more flexible alternative (Baugh 2006). These models use simple 
prescriptions to describe the various physical processes governing 
the evolution of the baryon content of the halos and aim to predict 
the fundamental properties of galaxies, such as their stellar mass 
and star formation history, ab initio. Adoption of a stellar popula¬ 
tion synthesis model and a choice of initial stellar mass function, 
allows one to translate these fundamental properties into directly 
observable properties. Although semi-analytic models still require 
observational data to constrain their parameters 1 , they have been 
shown to make realistic predictions for the evolution of the global 


1 Unlike HODs, the observational data used to constrain semi-analytic 
models need not be from the particular epoch of interest. Semi-analytic 
models can be constrained using local observations and then be used to 
make high redshift predictions. 


galaxy population out to high redshift (e.g. Lacey et al. 2011). How¬ 
ever, given the extreme cosmic evolution of galaxy clusters, which 
account for only a few per cent of all mass today, making accu¬ 
rate predictions for the properties of high redshift cluster galaxies 
remains a challenge as we shall see. 

Here we examine the predictions made by a semi-analytical 
galaxy formation model for the near-infrared (near-IR) photometry 
of galaxies in clusters at redshifts z > 1. Our decision to examine 
the statistics of clusters at this epoch in the near-infrared is mo¬ 
tivated by the aims of the Euclid mission, which will provide a 
deep survey over 15,000 deg 2 of the sky to a photometric depth 
of H < 24. As such, Euclid is predicted to provide a uniformly 
selected sample of approximately 60,000 clusters with a signal-to- 
noise greater than 3, with approximately 10,000 of these lying at 
z > 1 (Laureijs et al. 2011). Examining the near-IR predictions 
for this epoch is therefore extremely timely in the preparation for 
Euclid. The statistics that we consider are the colour-magnitude re¬ 
lation (CMR), due to its important role in cluster identification, and 
the cluster galaxy luminosity function (CGLF), which is one of the 
simplest statistics that can be made for the population of cluster 
galaxies. 

In §2 we introduce the galaxy formation model and describe 
how we select galaxies in clusters. The set of observed clusters, 
which we compare with the model predictions, are introduced in 
§3. In §4 we compare the model predictions for the CGLF, to ob¬ 
servational estimates between redshifts z = 1.2 and z = 1.6 and 
examine possible factors that might be causing the discrepancy be¬ 
tween the observations and the model predictions. Next, in §5, we 
compare the model predictions for the CMR, of cluster galaxies 
with the observed one. In §6 we examine the effect of varying se¬ 
lected model parameters. Finally, we draw our conclusions in §7. 

All synthetic and observed magnitudes have been converted to 
the AB system. 


2 GALAXY FORMATION MODEL 

In this section we describe the galaxy formation model that we em¬ 
ploy, starting with the N-body simulation used (§2.1) and followed 
by the semi-analytical model (§2.2). We then discuss the dust at¬ 
tenuation treatment used in the model (§2.3) and explain how we 
define a galaxy cluster (§2.4). 

2.1 N-body simulation 

The cosmological simulation that we use is a revision of the Mil¬ 
lennium Simulation (Springel et al. 2005), constructed using us¬ 
ing a cosmology consistent with the 7 year results of the Wilkin¬ 
son Microwave Anisotropy Probe (WMAP7, Komatsu et al. 2011). 
The cosmological parameters are: a baryon matter density Db = 
0.0455, a total matter density — Hb + Ocdm = 0.272, 
a dark energy density £2 a = 0.728, a Hubble constant Ho — 
100/ikms -1 Mpc -1 where h — 0.704, a primordial scalar spec¬ 
tral index n s = 0.967 and a fluctuation amplitude ag = 0.810. We 
shall refer to this simulation as the MS-W7 Simulation (Guo et al. 
2013). 

The hierarchical growth of cold dark matter structure is fol¬ 
lowed from redshift z — 127 to the present day, in a cubic volume 
of size 500/i _1 Mpc on a side. Halo merger trees are constructed 
using particle and halo data stored at 62 fixed epoch snapshots, 
which are spaced approximately logarithmically in expansion fac¬ 
tor. Details regarding construction of the halo merger trees can be 
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found in Merson et al. (2013) and Jiang et al. (2014). The MS-W7 
simulation uses 2160 3 particles to represent the matter distribution, 
with the requirement that a halo consists of at least 20 particles for 
it to be resolved. This corresponds to a halo mass resolution of 
Mhaio,iim = 1-87 x lO lo /i _1 M 0 , significantly smaller than ex¬ 
pected for the Milky Way’s dark matter halo. 

2.2 The GALFORM semi-analytical model 

To model the star formation and merger history of galaxies we 
use the GALFORM semi-analytical model of galaxy formation (Cole 
et al. 2000). The model populates dark matter halos with galaxies 
by using a set of coupled differential equations to determine how 
the various baryonic components of galaxies evolve (Baugh 2006; 
Benson 2010). 

For the work presented here we use a development of 
GALFORM that accounts for the following physical processes: (i) 
the collapse and merging of dark matter (DM) halos, (ii) the shock¬ 
heating and radiative cooling of gas inside DM halos, leading to the 
formation of galactic discs (iii) quiescent star formation in galactic 
discs, explicitly following the atomic and molecular gas compo¬ 
nents (Lagos et al. 2011, 2012), (iv) feedback as a result of su¬ 
pernovae, active galactic nuclei (Bower et al. 2006) and photo¬ 
ionisation of the inter-galactic medium, (v) chemical enrichment 
of stars and gas, (vi) dynamical friction driven mergers of galax¬ 
ies within DM halos, capable of forming spheroids and triggering 
starburst events, and (vii) disk instabilities, which can also trigger 
starburst events. 

Most of the published versions of GALFORM adopt a single 
Kennicutt (1983) initial mass function (IMF, see Baugh et al. 2005 
for an illustration of using a top-heavy IMF in starbursts) and up¬ 
dated versions of stellar population synthesis (SPS) models from G. 
& Chariot (1993). (See Gonzalez-Perez et al. 2014 for a compari¬ 
son of coupling GALFORM with different SPS models). By com¬ 
bining the star formation histories of the galaxies with the SPS 
models, GALFORM is able to calculate spectral energy distributions 
(SEDs) for the galaxies. Absolute magnitudes in a given photomet¬ 
ric band can be obtained by integrating the SED with the corre¬ 
sponding frequency-dependent filter response curve. (To calculate 
magnitudes in the observer-frame, a frequency shift is first applied 
to the filter response curve). All magnitudes and colours are total 
magnitudes. Apparent magnitudes are calculated using the redshift 
of the simulation snapshot to determine the distance modulus (see 
Merson et al. 2013). Note that the model magnitudes do not include 
any photometric uncertainties (see Ascaso et al. 2015 for a discus¬ 
sion of the impact of photometric errors on the colour-magnitude 
relation). 

GALFORM is able to track the global metallicity for the stars, 
as well as the hot and cold gas in the galaxy. Chemical enrich¬ 
ment is modelled using the instantaneous recycling approximation, 
with an effective yield and a recycled fraction that depend upon the 
choice of IMF. The yield is modified accordingly by metal ejection 
and feedback and hence is a function of the depth of the potential 
well of the galaxy. The rate at which gas is ejected from the galaxy 
due to supernovae explosions, M e j ec t, is given by, 

Meject - ( —) h0t M*, (1) 

V^disc / 

where Udi sc is the circular velocity of the galaxy disc at the half 
mass radius, M* is the star formation rate and cthot and Uhot are 
free parameters that govern the strength of supernovae feedback. 

In GALFORM, feedback due to AGN is implemented in halos 


that are undergoing quasi-static cooling, where, at fixed radius, the 
cooling time of the hot halo gas, t coo i, exceeds the dynamical free- 
fall time of the gas, m . Therefore feedback due to AGN can only 
occur when the condition, 



is satisfied, where a coo i is a free parameter. This condition is eval¬ 
uated at the cooling radius, r CO oi, which is defined, for a halo of a 
given age, as the radius at which the hot gas has only just had suf¬ 
ficient time to cool and collapse onto the galactic disc. Reducing 
the value of ct CO oi raises the minimum halo mass at which a quasi¬ 
static halo is established, thus allowing star formation to continue 
for longer in more massive halos. 

The free parameters in the Gonzalez-Perez et al. model were 
calibrated to reproduce the bj and K-band luminosity functions at 
z = 0 as well as to predict a reasonable evolution for the rest- 
frame K-band and UV luminosity functions. We stress at this point 
that this model has not been explicitly constrained using any obser¬ 
vations of high redshift clusters. 

Overall, the GALFORM model is able to make predictions for 
numerous galaxy properties, including luminosities over a substan¬ 
tial wavelength range extending from the far-UV through to the 
sub-millimetre. However, matching precisely the observed colour 
distribution of galaxies in the local Universe remains difficult for 
semi-analytical models (e.g. Font et al. 2008; Gonzalez et al. 2009; 
Guo et al. 2011). 

2.3 Modelling dust attenuation 

The attenuation by dust of the starlight from galaxies is modelled 
in GALFORM using a physically motivated method based upon 
the results of the radiative transfer code of Ferrara et al. (1999). 
The method assumes the dust to be distributed in dense molecular 
clouds embedded in a diffuse component. In this model the dust at¬ 
tenuation varies self-consistently with other galaxy properties, such 
as size, gas mass and metallicity, which are predicted by GALFORM 
(see also Fontanot et al. 2009b; Fontanot & Somerville 2011). The 
V-band optical depth when looking face-on through the centre of a 
galaxy, Ty, is assumed to be, 

0 M-coidZcoid 

Ty OC - 2 -5 ( 3 ) 

rp* 

disc 

where M co ia is the mass of cold gas in the galaxy (both atomic and 
molecular), Z co ia is the metallicity of the cold gas content of the 
galaxy and raise is the radius of the galactic disc. Given an extinc¬ 
tion law, the Ferrara et al. (1999) model provides dust attenuation 
factors as a function of wavelength, galaxy inclination, the ratio 
of bulge to disc radial dust scale length, r e /hR , the ratio of dust 
to stellar vertical scale heights, h Zj dust / h z , stars, and Ty (see Cole 
et al. 2000 and Gonzalez-Perez et al. 2014 for further details). 

2.4 Cluster galaxy selection 

When comparing the semi-analytic predictions to the observational 
estimates, we consider clusters of galaxies to be hosted by dark 
matter halos with mass, Mhaio, 

Mhalo ^ 1.2 x lO 14 /* - ^©. (4) 

This value is chosen as a compromise to ensure that we have a 
sufficiently large sample of halos that are massive enough to ade¬ 
quately represent our set of observed clusters. After imposing the 
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halo mass limit we are left with 43 halos at z ~ 1.4, 98 halos at 
z ~ 1.2 and 10 halos at z ~ 1.6. Cluster member galaxies are 
then taken to be those galaxies sharing a common host dark matter 
halo. Additionally, we place an aperture of 120 // by first assuming 
an observer placed at infinity, viewing the halos along the Carte¬ 
sian Z-axis of the simulation box. We use the positions of galaxies 
labelled by GALFORM as being central galaxies, which are located 
at the centre of mass of the most massive sub-halo of a halo, as the 
positions of the halo centres. We then apply the aperture using the 
projected distances between the galaxies and the halo centre. The 
effect of our choices for the halo mass limit and aperture size, as 
well as our method for modelling dust attenuation, are discussed 
further in §4.2. 

We stress that no colour selection is applied when selecting 
the GALFORM cluster galaxies and that, since we know the halo 
membership of the galaxies, the model predictions do not include 
contamination from foreground or background interlopers. For this 
work we are therefore examining the properties of model galaxies 
that truly are cluster galaxies, i.e. are hosted by cluster-mass halos 
selected according to our threshold of Mhaio > 1.2 x 1O 14 /i - 1 M 0 . 
An assessment of the effect of member incompleteness and in¬ 
terlopers, where cluster galaxies are selected according to their 
colours, is left for future work. 


3 OBSERVATIONAL DATASETS 

Here we briefly introduce the observational datasets against which 
we will compare the model predictions. 

Observations of clusters at high redshift, z > 1, are still lim¬ 
ited to of the order of a few tens of clusters, often with masses 
exceeding 10 14 /i _ 1 M©. Here, we consider a compilation of clus¬ 
ters at redshifts between z ~ 1.2 and z ~ 1.6, for which multi¬ 
band photometry is available. In particular, we focus on clusters 
for which photometry is available such that they could be identified 
by the Euclid mission. The clusters that we consider are listed in 
Table 1. 

Typically, the colours used to examine cluster galaxies are 
chosen such that the pairs of photometric bands bracket the red- 
shifted 4000A break, the strength of which is typically used as a 
proxy for the age of a galaxy and as a way of distinguishing pas¬ 
sively evolving galaxies from those that are undergoing star forma¬ 
tion when combined with other indexes (Kauffmann et al. 2003; 
Kriek et al. 2006, 2011), as illustrated in Fig. 1 for two synthetic 
galaxies at z — 1.4. At z ~ 1.2 the 4000 A break is shifted 
to 8800A, which lies between the i-band and z-band, whereas at 
z ~ 1.4 the 4000 A break, now shifted to 9600A, is bracketed by 
the z-band and J-band. At z ~ 1.6 the 4000A break is shifted to 
10400 A and is still bracketed by the z-band and J-band. 

For the majority of the clusters, mass estimates are available. 
In Table 1 we provide M 200 and M500 mass estimates from Ettori 
et al. (2004) and Jee et al. (2011). In many cases the mass esti¬ 
mates from gravitational lensing are larger than the mass estimates 
based upon the X-ray emission from the cluster, which is likely to 
be due to projection effects in the lensing estimate. However, X-ray 
mass estimates are also uncertain due to the limitations in our un¬ 
derstanding of the conditions in X-ray gas (e.g. Angulo et al. 2012). 

In the remainder of this section we provide further details on 
each of the galaxy cluster datasets. 

• Mei et al (2006b): Mei et al. present 1775 (F775W) and 
Z 850 (F850LP) observations of the clusters RX J0849+4452 and 


3000 4000 5000 6000 7000 



Ajbs a t z = 1.4 (A) 


Figure 1 . Examples of intrinsic spectra for two synthetic galaxies at redshift 
z = 1.4, over-plotted with transmission profiles for the i, z, J and H bands, 
as shown by the shaded regions. No dust attenuation is considered. Both of 
the galaxies were assumed to undergo a single instantaneous burst of star 
formation, one 300Myr ago (thin blue line) and the other 3Gyr ago (thick 
red line), prior to z = 1.4. (Both galaxies are assumed to have a metallicity 
of Z = 0.008). The spectra were generated using the PEGASE . 2 code 
(Fioc & Rocca-Volmerange 1999), assuming a Kennicutt (1983) IMF. The 
dotted, vertical line indicates the rest-frame wavelength of 4000A. 


RX J0848+4453, which together make up the Lynx Superclus¬ 
ter, located at z ~ 1.2. Applying the colour selection 0.8 < 
(i 775 — Z 850 ) < 1.1, the flux selection 21 < zsso < 24 and 
an aperture selection of 120", left 40 galaxies, of which 14 are 
confirmed cluster members and 26 are cluster member candidates. 
Galaxy colours were measured within the effective radii of the 
galaxies. Mei et al. assumed reddening due to dust to be described 
by a Calzetti et al. (2000) extinction law of E(B — V) = 0.027 
with Ai 775 = 0.054 and A Z850 = 0.040. 

• Strazzullo et al. (2006): Strazzullo et al. present K s -band 
imaging for the clusters RDCS J0910+5422 (K s < 21.5) and 
RDCS J1252.9-2927 (K s < 24.5), as well as Hi 60 (F160W) imag¬ 
ing for RX J0848+4453 (Hi6o < 25). They present an estimate for 
the CGLF for each cluster, as well as a composite estimate for all 
three of the clusters. The apertures applied are listed in Table 1. 

• Hilton et al. (2009): Hilton et al. present J-and K s -band pho¬ 
tometry, along with z 8 so — J and zsso — K s colour information, for 
64 galaxies selected as members of the cluster XMMXCS J2215.9- 
1738. The radial distance of the galaxies relative to the cluster X- 
ray source position extends out to 922 kpc. Spectroscopic redshifts 
are obtained for 24 of these galaxies. The photometry was corrected 
for Galactic extinction using the dust emission maps of Schlegel 
etal. (1998). 

• Strazzullo et al. (2010): Strazzullo et al. present multi¬ 
wavelength data for the cluster XMMU J2235-2557, one of the 
most massive virialised structures found beyond z ~ 1. They 
present estimates for the CGLF of the cluster in z 85 o, Hi 6 o and 
K s , with 10 a completeness limits of Z 775 < 25.3, H 160 < 25 and 
K s < 23. They also apply an aperture and select only those galax¬ 
ies within 83" of the cluster centre. The photometry was corrected 
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Table 1 . High redshift observational galaxy cluster datasets used for comparison with model predictions. Columns show: (i) Reference and cluster ID, (ii) 
estimated cluster redshift, (iii) photometry available, (iv) aperture used to identify cluster members, (v) estimated mass of cluster from X-ray measurements, (vi) 
estimated mass of cluster from weak lensing measurements. Cluster mass estimates are provided where known, with values converted to units of 10 14 /i _ 1 Mq, 
using the values of h specified by the original authors. Note that cluster RX J0848+4453 is listed twice. 


Cluster 

Cluster 

Photometry 

Aperture 

X-ray Mass 

Lensing Mass 


Redshift 


(Arcseconds) 

( 10 14 /i 1 M!q) 

(lO^/l^M©) 

Mei et al. (2006b) 






RX J0849+4452 

1.26 

1775, z 850 

120 

2.0 d= 1.0 a 

3.0 8 ™ 

RX J0848+4453 

1.27 

1775, z 850 

120 

0.96 d= 0.69 a 

9 9 + 0.7 b 
— 0.6 

Strazzullo et al. (2006) 






RDCS J0910+5422 

1.106 

Hieo,K 

73 

2.1 db 1.3 a 

q cr+ 0.8 b 

O.O — 07 

RDCS J1252.9-2927 

1.237 

Hieo,K 

59 

1.11 d= 0.24 a 

4.8±°', 6 

RX J0848+4453 

1.273 

Hieo,K 

47 

0.96 d= 0.69 a 

9 9+0.7 b 
^ - 0.6 

Hilton et al. (2009) 






XMMXCS J2215.9-1738 

1.46 

Z850,J,K, 

153 c 

i ^+ 0.86 b 
x -^-0.42 

3.0±H b 

Strazzullo et al. (2010) 






XMMU J2235-2557 

1.39 

z 850, Hi60, J, K, 

83 

A q+ 1-0 b 
4 ' d - 0.8 

c i +1-2 b 

°- i -i.o 

Snyder et al. (2012) 






ISCS J1426.1+3403 

1.136 

Hi60, 1814 

122 

- 

- 

ISCS J1426.5+3339 

1.163 

Hi60, 1814 

122 

- 

- 

ISCS J1434.5+3427 

1.243 

1775, Hi60 

120 

- 

i o + l -6 b 
i ‘°- 0.8 

ISCS J1429.3+3437 

1.262 

z 850, Hi60 

120 

- 

3.8lJ;p 

ISCS J1432.6+3436 

1.349 

z 850, Hi60 

119 

- 

q 7 + 1-8 b 
o.t _i 2 

ISCS J1433.8+3325 

1.369 

z 850, Hi60 

119 

- 

- 

ISCS J1434.7+3519 

1.372 

z 850, Hi60 

119 

- 

9 n+ 2.1 b 
z - u — 1.0 

ISCS J1438.1+3414 

1.413 

z 850, Hi60 

119 

o 9+3.7 b 
— 1.0 

9 9 + 1-8 b 
— 1.0 

Fassbender et al. (2014) 
XDCP J0044.0-033 

1.58 

J, K, i, V 

30 

2.1 

- 


a M 500 estimate from Ettori et al. (2004) 
b M 2 oo estimate from Jee et al. (2011) 

c Based upon maximum radial distance quoted in Table 1 of Hilton et al. (2009) 


for Galactic extinction using the dust emission maps of Schlegel 
et al. (1998). Strazzullo et al. argue that for their dataset the popu¬ 
lation of brightest galaxies in the inner region of the cluster can be 
considered to be quite well established, with about seventy per cent 
of these galaxies having measured spectroscopic redshifts. 

• Snyder et al. (2012): Snyder et al. present Hubble Space 
Telescope follow up observations of clusters selected from the 
Spitzer/IRAC Shallow Cluster Survey (ISCS, Eisenhardt et al. 
2008). They provide photometry for these clusters in the Hi 6 o-hand 
and one of the Isi 4 , i 775 or zsso bands. They report that the detec¬ 
tion catalogues for each cluster are more than 90 per cent com¬ 
plete for Hi 60 < 23.5. To identify cluster galaxies lying on the red 
sequence, they subtract from the colour of each galaxy a fiducial 
evolved CMR model for the Coma cluster and apply a selection 
based upon the residual. Reddening due to dust is assumed to be 
small; approximately E(B — V) <0.06. 

• Fassbender et al. (2014): Fassbender et al. present 
YLT/HAWK-I J and K s -band observations of the cluster XDCP 
J0044.0-033, complemented with the V and i-bands from Subaru 
archival imaging. They report that their observations are 100 per 
cent complete down to J < 23.9 and K s < 23.8. To maximise 
the signal-to-noise ratio of the cluster members to the interlopers, 
Fassbender et al. apply an aperture of 30 // . They estimate that the 
galaxies within this radius are 90 per cent cluster-associated. 


4 THE CLUSTER GALAXY LUMINOSITY FUNCTION 

In this section we compare observational estimates for the CGLF, 
which describes the number of galaxies per cluster as a function 
of apparent magnitude, with the predictions from our reference 
model. We stress again that for the model predictions, knowledge 
of the halo membership of the galaxies means that we can sim¬ 
ply select cluster galaxies using a halo mass selection and that no 
further colour selections are imposed. The observational estimates 
and semi-analytical predictions for the z-band CGLF are shown in 
Fig. 2, the H-band CGLF in Fig. 3 and the K-band CGLF in Fig. 4. 
The observational estimates are shown by the various data points, 
whilst the model predictions are indicated by the solid lines (with 
shaded regions indicating the Poisson uncertainties). 

4.1 Observational estimates 

In Fig. 2, Fig. 3 and Fig. 4 we show observational estimates for 
the CGLF from Strazzullo et al. (2006), Strazzullo et al. (2010) 
Fassbender et al. (2014), as well as our estimates based upon the 
datasets of Mei et al. (2006b), Hilton et al. (2009) and Snyder et al. 
( 2012 ). 

Strazzullo et al. (2006) and Strazzullo et al. (2010) estimate 
the CGLF for their galaxy cluster sample by using a reference field, 
down to an equivalent photometric depth or deeper, to statistically 
remove the contribution from background galaxies, which is known 
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Figure 2. z-band cluster galaxy luminosity functions (CGLF) at redshifts z ~ 1.2 (left) and z ~ 1.4 (right). The blue, solid line shows the predicted CGLF 
for the reference model, with dust attenuation. The shaded region indicates the Poisson uncertainty on this prediction. The dashed line shows the prediction for 
the reference model with no dust attenuation applied, i.e. using the fluxes intrinsic to the galaxies. The dotted line shows the prediction for the reference model 
when a dust attenuation similar to a Calzetti et al. (2000) extinction law with E(B — V) = 0.02 is assumed. All model predictions correspond to the CGLF 
for all cluster galaxies within an aperture of 120", hosted by dark matter halos of a i 0 ^ 1.2 x 10 14 A -1 Mq. Observational estimates of the luminosity 
functions are indicated by various data points. 




Figure 3. H-band cluster galaxy luminosity functions (CGLF) at redshifts z ~ 1.2 (left) and z ~ 1.4 (right). All model predictions correspond to the CGLF 
for all cluster galaxies, with galaxies selected in the same way as in Fig.2. Observational estimates of the luminosity functions are indicated by various data 
points. 
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Figure 4. K-band cluster galaxy luminosity functions (CGLF) at redshifts z ~ 1.2 (left) and z ~ 1.4 (middle) and z ~ 1.6 (right). All model predictions 
correspond to the CGLF for all cluster galaxies, with galaxies selected in the same way as in Fig.2. Observational estimates of the luminosity functions are 
indicated by various data points. 


to bias CGLF estimates (Andreon et al. 2005), particularly at bright 
magnitudes where the statistics are typically quite poor. The CGLF 
is estimated by subtracting the counts in the reference field (nor¬ 
malised to the solid angle of the cluster) from the counts in the 
cluster field. The uncertainties on the CGLFs estimated by Straz- 
zullo et al. (2006) and Strazzullo et al. (2010) are Poisson, with 
both clusters and possible background field interlopers summed in 
quadrature. The limits on the excess counts are determined based 
upon the upper and lower limits on the number of sources with a 
spectroscopic redshift or a photometric redshift within 3cr of the 
cluster redshift. Strazzullo et al. (2006) argue that the background 
contamination from lensed galaxies magnified by the cluster itself 
is small. They conclude that their estimate for the CGLF is con¬ 
sistent with previous determinations at similar or lower redshifts. 
Fassbender et al. adopt the same selection procedure as Strazzullo 
et al. (2006) when estimating the z ~ 1.6 K-band CGLF. 

For the datasets of Mei et al., Hilton et al. and Snyder et al. we 
have made simple estimates for the CGLF by counting all of the 
galaxies regarded as being cluster members by the original authors. 
From the Mei et al. data we estimate the z-band CGLF at z ~ 1.2. 
From the Hilton et al. data we estimate the z and K-band CGLFs at 
z ~ 1.4. From the Snyder et al. data we estimate the z and H-band 
CGLFs at z ~ 1.2 and z ~ 1.4. We apply no further selection 
beyond those placed originally by the authors and provide simple 
Poisson uncertainties on these counts. As such, our estimated un¬ 
certainties on these observational luminosity functions may well be 
under-estimated and our CGLF estimates may be biased by mem¬ 
bership incompleteness. 

At z ~ 1.4, in the z and H bands, we compare our estimates 
for the CGLF to the estimates from Strazzullo et al. (2010), though 
these are for a slightly more massive cluster and adopt a smaller 
aperture when selecting the member galaxies (as we demonstrate 
in §4.2.1, changing the aperture size has a greater impact at the 
faintest magnitudes). We find that our CGLF estimates from the 
Snyder et al. data are consistent within error with the estimates 
from Strazzullo et al.. When constructing their dataset, Snyder et al. 
identified galaxies as cluster members based upon the likelihood of 
the galaxy lying on the cluster RS. For a subset of the galaxies these 
authors were able to confirm that they indeed are on the RS of the 


clusters. In Fig. 2 and Fig. 3 we plot the Snyder et al. estimate for 
the CGLF using only the confirmed RS galaxies and can see that 
this is in excellent agreement with the estimate using the full Sny¬ 
der et al. dataset, suggesting that our CGLF estimates are not being 
significantly biased by possible interloping non-member galaxies. 

At z ~ 1.2, there are no existing estimates of the CGLF in 
the H or z-bands against which to compare ours. In the these bands 
our estimates from the Snyder et al. dataset show a bright-end fall 
off similar in shape to the CGLF estimates at z ~ 1.4. At the 
faintest magnitudes, however, the CGLFs show a down-turn. Sny¬ 
der et al. report 5 a detection limits of 26.0 — 26.3 in the z-band and 
24.4 — 24.8 in the H-band, with 90% completeness at H ~ 23.5. 
It is possible therefore that the observed down-turns may be a re¬ 
sult of incompleteness. When determining their samples of cluster 
galaxies, Snyder et al. place a selection limit of H ~ 22 — 23, de¬ 
termined for each cluster by evolving the characteristic brightness 
of the Coma cluster to the redshift of each cluster. This selection 
may also be contributing to the down-turns. 

In the z-band, our estimate for the CGLF from the Mei et al. 
shows good agreement with the Snyder et al. estimate for magni¬ 
tudes bright-wards of z ~ 22.2. Faint-wards of this value, however, 
the estimates diverge. We note that these estimates are based upon 
galaxies from only one or two clusters, so discrepancies are indeed 
possible due to cosmic variance across the cluster populations. An¬ 
other possible cause could be due to the i — z colour selection that 
Mei et al. place in order to identify early-type galaxies. We note, 
however, that the uncertainties on our CGLF estimates are simple 
Poisson errors and so may well be underestimates of the true un¬ 
certainties. If this is the case, then it is possible that our estimates 
for these two datasets are consistent within error. 

Our estimates for the z and K-band CGLFs from the Hilton 
et al. dataset are generally consistent with the other estimates pre¬ 
sented, though there are one or two bins for which the counts are 
higher than the other estimates, for example at z ~ 23.2 and 
K ~ 21.2. Just under two-thirds of galaxies in the Hilton et al. 
data have photometric redshifts only and so this variation could be 
caused by interlopers. In addition, variation due to cosmic variance 
between CGLF measurements for different clusters would be ex¬ 
pected. 
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Figure 5. H-band CGLFs at z ~ 1.4 as predicted by GALFORM assum¬ 
ing different apertures for selecting the member galaxies of each halo, as 
indicated in the legend. For each choice of aperture a halo mass limit of 
A^haio ^ 1.2 x 10 14 /i _1 Mq was adopted. The data points are the same 
as in the right-hand panel of Fig. 3. 


4.2 Semi-analytical predictions 

The predictions for the CGLF from our reference GALFORM model 
are shown in Fig. 2, Fig. 3 and Fig. 4 by the solid blue lines, with 
shaded regions indicating the Poisson uncertainties. These predic¬ 
tions correspond to the CGLF for all cluster galaxies, i.e. all galax¬ 
ies in halos with mass Mhaio ^ 1.2 x 10 14 /i _1 M q , within a pro¬ 
jected aperture of r < 120". 

There is a large discrepancy between the model predictions 
and the observational estimates around the knee of the CGLF. This 
discrepancy is most obvious in the H and K-band CGLF compar¬ 
isons, though it is still apparent in the z-band CGLF, particularly 
at 2 ~ 1.4. In the brightest and faintest magnitude bins, however, 
the model predictions for the CGLF are broadly consistent with the 
observational estimates. 

In the H and K-bands, the model predictions for the CGLF 
show a faint-end slope that is quite flat, in agreement with the faint- 
end slope seen in the observations. In the K-band, however, the 
normalisation of the slope in the model predictions is lower than 
that of the observations, particularly at z ~ 1.2. 

We now consider three factors that could affect our compar¬ 
ison with cluster data, particularly our estimation of the CGLF. 
These are: the size of the aperture we apply to the halo (§4.2.1), 
our choice of halo mass limit applied to the model (§4.2.2) and our 
modelling of attenuation due to dust (§4.2.3). 


4.2.1 Aperture size 

We select the GALFORM cluster galaxies using an aperture size of 
r < 120 // , which is consistent with the apertures placed by Mei 
et al. (2006b) and Snyder et al. (2012). Given the cosmology used 
in the MS-W7 simulation, an aperture of 120" corresponds approx¬ 
imately to a distance of 500/i _1 kpc at z ~ 1.4. Varying the aper¬ 
ture size will change the number of satellite galaxies that are in¬ 


cluded. We must therefore examine how a change in aperture size 
affects our estimates of the CGLF. 

In Fig. 5 we plot the predictions for the reference model for the 
H-band CGLF at 2 ~ 1.4 when keeping the halo mass limit fixed 
and allowing the aperture size to vary. The variation of the aper¬ 
ture size induces a change in the abundance of cluster members, 
particularly for the faintest magnitude bins where a change in the 
aperture size will lead to different numbers of faint satellite galax¬ 
ies being selected. This change in the abundance, typically within 
a factor of two, is consistent with observational uncertainties. All 
of the apertures we consider correspond to distances smaller than 
the virial radius of a typical cluster-sized dark matter halo and so 
remove the most distant satellites from the comparison such that 
we are comparing only the cores of the clusters. At brighter magni¬ 
tudes, H < 21.5, the change in the abundance becomes smaller due 
to the increasing number of bright, central galaxies, which will al¬ 
ways be selected. If we remove the aperture altogether then we see 
an increase of approximately a factor of two in the abundance of 
cluster members with H > 21.5, as all of the most distant satellites 
are now included (shown by the red line in Fig. 5). 

From Fig. 5 we can see that an aperture of 120" is consistent 
with many of the aperture choices in Table 1 and provides a suitable 
match to the counts just faint-wards of the break in the CGLF. 

4.2.2 Halo mass limits 

We select cluster galaxies in GALFORM as galaxies hosted by halos 
above a threshold mass of 1.2 x 10 14 /i _1 M q . Using a fixed aper¬ 
ture of 120" we find that variation of the halo mass limit between 
1.0 x 10 14 /i -1 M© and 1.8 x 10 14 /i _1 M© produces a negligible 
change in predicted H-band CGLF at z ~ 1.4. (Examination of 
larger halo masses is not possible due to the limited volume of the 
MS-W7 simulation.) We find that the choice of halo mass limit has 
a negligible effect on the CGLF. 

4.2.3 Attenuation due to dust 

The dust content of a galaxy can have a drastic effect upon the ob¬ 
served colour of the galaxy. In Figs. 2, 3 and 4 we compare the 
model CGLF, both with and without dust attenuation, with the ob¬ 
servations. In the reference model dust attenuation has a large im¬ 
pact around the knee of the CGLF. At the knee of the CGLF the 
reference model is about a factor of 4.5 below the observations. 
The observed CGLF can be reproduced if instead we consider a 
model without dust attenuation. However, such a model predicts an 
unrealistic luminosity function at z = 0, as we shall discuss later. 

To help understand this result, we compare the intrinsic and 
attenuated magnitudes of the GALFORM galaxies in our reference 
model. We define the difference, 

AM = M a tt. - Mint., (5) 

where Mint is the intrinsic, dust-free magnitude of a GALFORM 
galaxy and M att . is the magnitude of this galaxy attenuated using 
the dust model described in §2.3. In Fig. 6 we plot the difference 
between these two magnitudes as a function of H-band intrinsic 
magnitude. The majority of galaxies display very little dust attenu¬ 
ation and so have a negligible difference between their intrinsic and 
attenuated magnitudes. 

There are, however, a small number of galaxies, with intrinsic 
magnitudes between 21 < Hint. < 23, that display dust attenu¬ 
ation larger than one magnitude. If we examine again the H-band 
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Table 2. Fits to the red sequences predicted by the Gonzalez-Perez et al. 
(2014) GAL FORM model for various colour-magnitude spaces between red- 
shifts z— 1.2 and z = 1.6 using subsets of model cluster galaxies selected 
by a cut in specific star formation rate, sSFR (a passively evolving sample). 
Note that the samples were additionally limited to galaxies brighter than 
25 th magnitude. 


Magnitude 

Colour 

Slope 

Zero-point 

REDSHIFT: 1.2 




z 

i — 

z 

-0.0118 ±0.0005 

1.0244 ±0.0009 

H 

i — 

H 

-0.146 ±0.002 

2.436 ± 0.003 

K 

J- 

K 

-0.118 ±0.001 

0.740 ± 0.002 

REDSHIFT: 1.4 




J 

z — 

J 

-0.060 ± 0.003 

1.348 ±0.004 

H 

z — 

H 

-0.115 ±0.004 

1.839 ±0.005 

K 

J - 

K 

-0.112 ±0.003 

0.849 ± 0.004 

REDSHIFT: 1.6 




K 

J - 

K 

-0.106 ±0.007 

0.99 ±0.01 


Figure 6. Difference in GALFORM H-band attenuated magnitudes, H a tt, 
and H-band intrinsic magnitudes, Hi nt , as a function of intrinsic magni¬ 
tude. (The default dust attenuation calculation is assumed). The colour bar 
indicates the number of galaxies in each pixel. The redshift, as well as the 
halo mass and aperture used to select the GALFORM galaxies, are shown in 
top right of the panel. 


tween the model and observations could, for example, impact upon 
the calibration of cluster-finding algorithms for next generation sur¬ 
veys such as Euclid. 


CGLF at z ^ 1.4 (right-hand panel, of Fig. 3), we can see that 
this magnitude range corresponds approximately to the knee of the 
CGLF, where the maximum discrepancy between the model and 
the observations occurs. For magnitudes outside this range, the at¬ 
tenuation is minimal. Hence, the size of the attenuation applied in 
our reference model is not constant with magnitude, unlike, for ex¬ 
ample, the simple Calzetti et al. (2000) fitting formula that is often 
applied to low redshift star-forming galaxies. With our reference 
attenuation calculation, the large dust attenuation in the brightest 
galaxies causes these galaxies to be shifted out of the brightest 
magnitude bins and to pile up in the magnitude bins faint-wards 
of the characteristic magnitude. Although these galaxies constitute 
only of the order 5 — 10 per cent of galaxies with H at t. < 25 for 
example, the rapidly declining number of galaxies in the brightest 
bins means that this magnitude shift has a significant impact upon 
the CGLF. 

As a final demonstration that these highly dust attenuated 
galaxies are the cause of the discrepancy between the model CGLF 
and the observations, we calculate again the CGLF for our refer¬ 
ence model but now assume a weaker dust attenuation, similar to a 
Calzetti et al. (2000) extinction law with E(B — V) = 0.02. The re¬ 
sulting CGLF predictions are shown as blue dotted lines in Figs. 2, 
3 and 4. In each instance, we see that weakening the attenuation in 
this way leads to a much better agreement between the model and 
the observations. 


5 THE CLUSTER COLOUR-MAGNITUDE RELATION 

Having examined the abundance of semi-analytical cluster galax¬ 
ies, we now consider their colours. Specifically, we consider the 
colour-magnitude relation (CMR). In the CMR the red sequence 
(RS) is often used observationally to detect clusters as it is expected 
to be dominated by early-type cluster members. It is desirable to 
examine the model predictions for the RS as any discrepancy be- 


5.1 Red sequence fitting 

To describe the location of the model RS, we use linear regression 2 
to fit the optimised slope, s, and zero-point, C22.5, for the relation, 

c = s( m - 22.5) + C 22 . 5 , (6) 

where m is the galaxy magnitude, e.g. H, and c is the corresponding 
colour, e.g. z — H. Since the RS is expected to be dominated by pas¬ 
sively evolving galaxies and since GALFORM is able to provide us 
with values for the star formation rates of galaxies, as well as which 
halos they belong to, we are able to fit a RS using just those galaxies 
that are passively evolving and reside in halos above our specified 
halo mass threshold. Therefore, we stress that this is not meant to 
mimic observational methods for determining the RS, but simply 
to provide an estimate of the RS predicted by the model. To deter¬ 
mine the passively evolving cluster galaxies, we apply a cut in in¬ 
stantaneous specific star formation rate (sSFR) and fit to only those 
galaxies that satisfy log 10 (sSFR/Gyr -1 ) ^ — 1, which provides a 
reasonable distinction between actively star-forming and passively 
evolving galaxies in observational data (e.g. Williams et al. 2009) 
and hydrodynamical simulations (e.g. Romeo et al. 2008; Furlong 
et al. 2015). We therefore regard the fit to the passively evolving 
galaxies as our measure of the true RS predicted by the model, i.e. 
the RS of those cluster galaxies that are truly ‘red and dead’. Note 
that in addition to the sSFR cut, we also select only those galaxies 
brighter than 25 th magnitude in the appropriate band (z, J, H or K). 

The fits to the passively evolved cluster galaxies provide an 
excellent description of the model RS down to 25 th magnitude, 
as evident from Figs. 7 and 8 where the fits to the passive RS are 
shown by the red, dashed lines. In these figures the greyscale pix¬ 
els show the distribution of all cluster galaxies in the model, pas¬ 
sive and star-forming, normalised by the number of clusters (i.e. the 

2 To determine the optimised value for the slope and zero-point we provide 
the list of galaxy magnitudes and colours to the curve_f it function in the 
Scientific Python library, scipy (http : //www. scipy. org/). 















10 Merson et al. 




o.io 

0.09 

0.08 

0.07 

0.06 

0.05 

0.04 

0.03 

0.02 


2.0 



0.5 


0.0 


-- log 10 (sSFR/Gyr ) ^ —1.0 
☆ XMMU J2235-2557 
# XMMXCS J2215.9-1738 (spec-z) * 
^ XMMXCS J2215.9-1738 (phot-z) 

iiiiiiiiiiiiiiiiiiiiiiiiiii.iii. i 


3.0 

2.5 

2.0 

K 1.5 

I 

N 

1.0 

0.5 

0.0 


M halfl > 1.2 x lO 14 ft _1 M 0 r < 120' 


1.4 



20.0 20.5 21.0 21.5 22.0 22.5 23.0 
J 


-- logio(sSFR/Gyr -1 ) ^ —1.0 * 

☆ XMMU J2235-2557 
| ISCSJ1433.8+3325 
i ISCSJ1432.6+3436 
| ISCS J1434.7+3519 
f ISCS J1438.1+3414 


20.0 20.5 21.0 21.5 22.0 22.5 23.0 
H 


Figure 7. GALFORM colour-magnitude relations at redshifts z ~ 1.2 (upper panels) and z ~ 1.4 (lower panels). The observed cluster galaxies are shown by 
the various data points. Greyscale pixels show the number of cluster galaxies per cluster as predicted by the fiducial GALFORM model. The red dashed lines 
show linear fits to the red sequence predicted by GALFORM. The halo mass lower limit and aperture, r, used to select the GALFORM galaxies are shown in the 
label at the top of each panel. 


number of halos above the halo mass threshold). In each case the 
distribution is dominated by a clear and well-defined RS, especially 
at faint magnitudes, with very little indication of a blue cloud. 


At bright magnitudes the model RS appears to display a 
prominent plume of galaxies with colours extending redwards 
above the RS. This plume is visible in many of the colour-spaces 
that we consider, in particular i — H, z — H and J — K. As we shall 
see in §5.3, the galaxies in the plume are star-forming and so are 
not identified by our sSFR selection. As such, they are not included 
in the fitting and do not bias the fits to the slope or the zero-point. 
In contrast, simply fitting the RS to a straight-forward flux-selected 
sample, i.e. without selecting just the passively evolving galaxies, 
would lead to biased fits. The results of the fits, i.e. the slope and 
intercept, as well as their corresponding uncertainties, are provided 
in Table 2. 


5.2 Comparison with observed clusters 

We now compare the GALFORM prediction for the cluster CMR 
with observational measurements at z > 1. In Figs. 7 and 8 we 
show the colours of the observed cluster galaxies on top of the 
GALFORM predictions, shown by the greyscale pixels. In addition, 
we show as a red, dashed line the fit to the passively evolving 
RS. The distribution of GALFORM cluster galaxies appears to be 
qualitatively in good agreement with the observations for each of 
the colour-magnitude spaces that we consider. In each instance the 
model RS is qualitatively in good agreement with observations. For 
all cases in Figs. 7 and 8, the scatter in the distribution of colours 
from GALFORM is consistent with or smaller than the spread in the 
observations, though we must recall that no photometric errors are 
included in the model predictions. 

Comparing our estimates for the slope and zero-point (see 
Eq. 6) of the model RS with fits available for the observed clus¬ 
ters, we find a reasonable agreement for the longer-wavelength 
colours, particularly for the zero-point estimates. For example, our 
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fit to the RS in z — J are consistent with the fit from Hilton et al. 
(2009) who estimated a slope of —0.049 zb 0.062 and a zero- 
point of 1.335 =b 0.046. For the clusters ISCS J1433.8+3325, ISCS 
J1432.6+3436, ISCS J1434.7+3519 and ISCS J1438.1+3414 Sny¬ 
der et al. determine a range of values for the z — H zero-point. 
Taking the mean of these values and adding the uncertainties in 
quadrature gives an estimate of 1.78 zb 0.09, which is consistent 
within la with our z — H fit. For i — z, our fit to the RS zero-point 
are close to the fit of Mei et al. (2006b), who found a zero-point of 
0.99 zb 0.01. Our fit to the slope, however, are shallower than their 
fit of-0.031 zb 0.012. 

5.3 Trends in galaxy properties 

We conclude our analysis of the cluster RS by examining the prop¬ 
erties of the galaxies in the CMR as predicted by our reference 
GALFORM model. In Fig. 9 we show a selection of the studied prop¬ 
erties for J — K. In these plots the colour-magnitude space has been 
divided into pixels and the colour maps show the median value of a 
particular property for all of the cluster galaxies that fall in that par¬ 
ticular pixel. We can see immediately a clear trend in galaxy prop¬ 
erties along the RS. The upper panel of Fig. 9 shows that the bright¬ 
est and reddest galaxies typically having higher SFRs than faint red 
galaxies. The middle panel also shows that the brightest and reddest 
galaxies are also typically the most metal rich. Given also that these 
bright, red galaxies also have the largest stellar masses, this result 
suggests a positive correlation between stellar mass, metallicity and 
SFR. The existence of such a correlation between models and ob¬ 
servations has been debated in the literature (e.g. Mannucci et al. 
2010; Magrini et al. 2012; Lilly et al. 2013; Obreja et al. 2014). 

In addition, the plume of bright, red galaxies that we have pre¬ 
viously commented on is clearly evident in Fig. 9. The galaxies in 
the plume are revealed to typically be highly star forming (with star 
formation rates in excess of 1 — 10 /i _1 M 0 yr _1 ), have reservoirs 
of cold gas much larger than other galaxies on the RS, be heav¬ 
ily attenuated by dust and have metallicities that are richer than 
the other galaxies on the RS. The lower panel of Fig. 9 shows the 
galaxies in the plume have very large optical depths, much larger 
than other cluster galaxies on the RS. The large dust attenuation 
for the galaxies in the plume would suggest that these are the same 
galaxies that are responsible for the discrepancy between the ob¬ 
served CGLF and the model prediction, though the exact cause of 
their large attenuation is not immediately clear. If we examine the 
predicted RS when we apply a weak Calzetti-like attenuation to the 
intrinsic magnitudes of the galaxies, then we find that the plume is 
removed. This happens at the expense of making the majority of 
bright galaxies up to one magnitude bluer, leading to a bluewards 
break in the bright RS. This suggests that the plume is an artefact 
of the model and that simply reducing the strength of the dust at¬ 
tenuation in the model is not an acceptable solution. Instead, this 
perhaps hints at an additional underlying problem in the model. 

Although we suspect that the plume of star-forming galaxies 
in the model RS is an artefact, the scatter in the observations hints at 
the existence of some observed cluster galaxies with colours as red 
as those in the plume. Several observations of our sample of high 
redshift clusters also suggest ongoing star formation activity. For 
the cluster XDCP J0044.0-2033, Fassbender et al. (2014) were able 
to obtain spectroscopic redshifts for a few galaxies with significant 
[Oil] emission, which is often taken as an indicator for ongoing 
star formation. One or two of these galaxies, which we have high¬ 
lighted in the bottom panel of Fig. 8, have very red J — K colours 
suggesting such red galaxies might indeed be found in clusters in 
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Figure 8. J — K vs. K colour-magnitude relations for cluster galaxies at 
redshifts z ~ 1.2 (top panel), 2 : ~ 1.4 (middle panel) and 2 : ~ 1.6 (bot¬ 
tom panel). The colours of the observed galaxies are shown by the various 
data points. The observations at 2 : ~ 1.6, from Fassbender et al. (2014), 
have been split according to distance r from the estimated cluster centre: 
galaxies within r < 13" (which typically have a 75% membership prob¬ 
ability) and galaxies with 13" < r < 30" (which typically have a 50% 
membership probability). Galaxies that have been spectroscopically con¬ 
firmed are highlighted as being passive (no Oil detection) or star-forming 
(Oil detection). Note that three of the spectroscopically confirmed galaxies 
lie at r > 30". Cluster data from Hilton et al. (2009) is split into galaxies 
with spectroscopic redshifts (spec-z) and those without spectroscopic red¬ 
shifts (photo-z). As before, the greyscale pixels show the prediction for the 
fiducial GALFORM model and the red, dashed line shows the linear fit to 
the predicted red sequence. The halo mass and aperture used to select the 
GALFORM galaxies are shown at the top of each panel. 
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Figure 9. Distribution of various galaxy properties across the J — K vs. 
K colour-magnitude space at z ~ 1.4. The panels, which show the dis¬ 
tribution of star-formation rate, stellar metallicity and combined half-mass 
radius of the disc and bulge of the galaxies, are labelled accordingly. Each 
pixel is coloured according to the median value for the galaxies that lie in 
that pixel, with the colour scale shown at the right of each panel. The halo 
mass and aperture used to select the galaxies are shown at the top of each 
panel. 



reality, though the large spread in the observations makes this un¬ 
clear. Demarco et al. (2007) presented spectroscopy for the cluster 
of Strazzullo et al. (2006) and measured [Oil] emission lines in 
38 cluster members. They estimated that the SFRs of those galax¬ 
ies are in the range 0.5 — 2 MQyr -1 , with the median SFR being 
« 0.7 M©yr _1 . Similarly, Strazzullo et al. (2010) reported SFRs 
of a similar magnitude for cluster galaxies from photometry, rest- 
frame FUV. All of these SFR tracers suggest inferred SFRs that are 
not as high as the SFRs predicted by the model for those galax¬ 
ies in the plume above the RS, although these SFR tracers can be 
heavily obscured. However measurements of near-IR spectroscopy 
targeting the Ha line, which is a more reliable SFR tracer, and IR 
photometry from Spitzer and Herschel hint to a higher SFR. For 
example, Valentino et al. (2015) find the H a luminosity of clus¬ 
ter CL J1449+0856, at 2 = 1.99, to be significantly higher than 
measurements in the field at the same epoch. They attribute this 
to an enhanced specific star formation rate in the cluster. In addi¬ 
tion, based upon measurements of the H a emission from galaxies 
in the Snyder et al. (2012) clusters, Zeimann et al. (2013) infer un¬ 
obscured SFRs of up to 200 M©yr -1 for galaxies right down in the 
cluster cores. SFR measurements of the same clusters from Spitzer 
24/xm observations, as well as Herschel SPIRE data, find similarly 
high SFRs, with the SFR in clusters increasing rapidly with redshift 
(Brodwin et al. 2013; Alberts et al. 2014). 


6 DISCUSSION 

We have seen in §4.2.3 that the dust attenuation calculation used in 
our reference GALFORM model predicts a large dust attenuation for 
cluster galaxies that leads to the CGLF predicted by the model be¬ 
ing inconsistent with observations. In addition, in §5 we have seen 
that this large attenuation also produces a plume of very red galax¬ 
ies above the RS, which could potentially bias predictions for the 
RS. We have seen that whilst applying a weaker dust attenuation re¬ 
moves the discrepancy between the observed CGLF and the model 
predictions, this leads to many of the galaxies in the RS being made 
too blue. Here we examine the cause of the large dust attenuation 
that produces the tension between the model predictions and obser¬ 
vations. 


6.1 Galaxy stellar mass-size relation 

From Eq. 3 we can see that the dust attenuation predicted in our 
reference model is affected by the predicted galaxy sizes, cold gas 
masses and the assumed distribution of dust with respect to stars. 
We have tested that changing the distribution of dust with respect 
to that of the galaxy stars has little impact on the CGLF. In Fig. 10 
we show how selected galaxy properties change with as a function 
of the difference between the attenuated and intrinsic H-band mag¬ 
nitudes of the GALFORM galaxies. We see from the bottom panel 
that galaxies with a large magnitude difference have a larger opti¬ 
cal depth We can also see that galaxies with large attenuation typi¬ 
cally have larger reservoirs of cold gas and smaller radii. Both the 
predicted galaxy sizes and cold gas masses are fundamental pre¬ 
dictions from the model that are directly related to the modelling 
of the cooling of gas and feedback processes (see Cole et al. 2000 
for details of how galaxy sizes are calculated in the model). As 
such, directly modifying these properties is a complex procedure 
and beyond the scope of this paper. We do, however, in §6.2 briefly 
explore how varying selected parameters of the reference model 
affects the predictions for the CGLF. 
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Figure 10. Difference between the H-band attenuated magnitudes and in¬ 
trinsic magnitudes, as in Fig. 6, but now showing the correlation with se¬ 
lected galaxy properties: cold (HI, HII and He) gas mass, half-mass radius 
and V-band optical depth (through centre of galaxy when face-on). The 
pixels are coloured according to the median value of the galaxy property 
for the galaxies in that pixel as shown by the key on the right side of each 
panel. The redshift, as well as the halo mass and aperture used to select the 
GALFORM galaxies, are shown in each panel. 
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Figure 11. Comparison of the mass-size relation predicted by GALFORM, 
with observational measurements for high-redshift early-type galaxies 
made by Papovich et al. (2012); Szomoru et al. (2012) and Saracco et al. 
(2014). The observational estimates were converted from effective radius, 
r e , to half-mass radius, rhaif> using the conversion, r*h a if = 1.35r e . The 
upper panel shows the model predictions and observational measurements 
at 2 : ~ 1.2, whilst the lower panel shows the predictions and measurements 
at 2 : ~ 1.6. The halo mass and aperture used to select the semi-analytical 
cluster galaxies are shown in the top -left-hand corner of each panel. 


Perhaps the easiest of these properties to compare against ob¬ 
servations is the size of the galaxies. In Fig. 11 we compare the 
stellar mass-size relation for the model cluster galaxies with ob¬ 
servational estimates from Szomoru et al. (2012) at 2 : ~ 1.2 and 
from Papovich et al. (2012) and Saracco et al. (2014) at z ~ 1.6. 
Note that the galaxies we show from Papovich et al. are those se¬ 
lected as cluster galaxies, whereas the galaxies from Szomoru et al. 
and Saracco et al. are from global samples 3 . For stellar masses 
M * < lO lo /i - 1 M 0 , the mass-size relation for the model, which is 
shown by the greyscale pixels, is approximately flat but with a large 


Some observations at high redshift appear to be consistent 
with high redshift galaxies having a small dust attenuation (e.g. 
Meyers et al. 2012). If we assume therefore that there is indeed 
negligible dust attenuation in high redshift cluster galaxies, then 
from Eq. 3 we can see that too large an optical depth might indicate 
that the predicted sizes of such galaxies are too small, the amount 
of cold gas in the galaxies is too high or the metallicity of the cold 
gas is too high, or a combination of all three. 


3 To convert the stellar mass estimates from Papovich et al. and Saracco 
et al. from the Chabrier initial mass function (IMF) to the Kennicutt (1983) 
IMF we adopt a conversion factor of -0.09 dex (Mitchell et al. 2013). To 
convert the stellar masses estimates of Szomoru et al., we first converted 
from the Kroupa (2001) IMF to the Salpeter (1955) IMF using a conversion 
factor of -1.6 dex (Fontana et al. 2004) and then converted from the Salpeter 
(1955) IMF to the Kennicutt (1983) IMF using a conversion factor of +1.4 
dex (Fontana et al. 2004). 
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scatter of typically 0.4 — 0.5 dex. The mass-size relation drops off 
towards the highest masses. Therefore, although there is agreement 
between the observations and the model predictions, the most mas¬ 
sive cluster galaxies, with stellar masses M* > 10 lo /i -1 M©, are 
typically smaller than the observed sizes, with some model galax¬ 
ies being up to an order of magnitude smaller. Examining how the 
stellar mass-size relation correlates with other galaxy properties, 
we can see in Fig. 12 that the population of very massive, com¬ 
pact galaxies in the model typically have the highest optical depth, 
which is not surprising given Eq. 3. We also find that these high- 
mass, compact galaxies have extremely large cold gas masses of 
M co id,gas ~ lO lo /i _1 M 0 (comparable to their stellar masses) as 
well as the highest star formation rates (between l/i _1 M©yr _1 
and lO/i -1 M 0 yr -1 ). 

Overall, the large reservoirs of cold gas and the compact sizes 
of the most massive cluster galaxies appear to be the cause of the 
large dust attenuation, which leads to the discrepancy between the 
observed CGLF and that predicted by the model. Our results there¬ 
fore suggest that the problem of the large dust attenuation is ac¬ 
tually due to an underlying problem of either the model under¬ 
predicting the sizes of the most massive cluster galaxies or the 
model allowing too much cooling of gas in these halos. Gonzalez 
et al. (2009) come to a similar conclusion when comparing pre¬ 
dictions from the Baugh et al. (2005) and Bower et al. (2006) 
GALFORM models with the colours of galaxies in the Sloan Digital 
Sky Survey (SDSS, York et al. 2000). Gonzalez et al. find that the 
model over-predicts the number of bright, blue galaxies and also 
that the bulge-dominated bright galaxies have sizes up to a factor 
of ten times smaller than SDSS galaxies of equivalent luminosity. 
Following an examination of several of the model parameters, they 
attribute the problem as being due to an over-simplified treatment 
of the sizes of galaxy merger remnants. At lower redshift, Weinzirl 
et al. (2014) compare the predictions of the semi-analytical model 
of Neistein & Weinmann (2010) with Hubble Space Telescope 
observations of the Coma Cluster and find that the model over- 
predicts the mass fraction of cold gas for galaxies in halos with 
Coma-like properties. 


6.2 Robustness of model predictions to parameter changes 

In the previous section we concluded that the discrepancy be¬ 
tween the observed CGFF and the model predictions is likely being 
caused by the model predicting too much cold gas and too small 
sizes for the most massive cluster galaxies. Since these two funda¬ 
mental galaxy properties are sensitive to many other model parame¬ 
ters, we now vary some of the key parameters that we would expect 
to have the greatest effect upon the high redshift cluster galaxy pop¬ 
ulation in order to gain some insight into the physics shaping the 
model predictions for cluster galaxies. We note that most of the 
parameters we choose to vary will have a greater impact on the 
cold gas masses of the galaxies rather than the galaxy sizes. For 
this exercise we will vary each parameter independently and use 
the prediction for the z ~ 1.4 CGFF as an indicator of possible 
improvements, since the discrepancy is most noticeable in the pre¬ 
dictions for the CGFF. We note, however, that for more a extensive 
search varying multiple parameters simultaneously and examina¬ 
tion of multiple galaxy statistics would be necessary. We leave such 
a search for future work. 

The main parameters that we expect to have the greatest influ¬ 
ence on the predictions for the cluster galaxy population are those 
governing the heating and cooling of gas in the most massive halos 
as well as those parameters governing treatment of galaxy merg- 
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Figure 12. Distribution of selected galaxy properties along the mass-size 
relation as predicted by GALFORM. The various properties: star formation 
rate, half-mass radius and V-band optical depth, are labelled in the corre¬ 
sponding panel. Note that the half-mass radius includes both the disc and 
bulge of the galaxy and the optical depth is the depth measured through 
the centre of the galaxy when face-on. The pixels are coloured according to 
the median value of the galaxy property for the galaxies in that pixel. The 
redshift, as well as the halo mass and aperture used to select the GALFORM 
galaxies, are shown in each panel. 


ers. Besides these parameters, we also examined varying the scale 
height of the dust in the model galaxies. Although this parameter 
cannot improve the model predictions for the galaxy sizes or cold 
gas masses, it may allow an improved recovery of the CGFF. We 
find, however, that this has limited impact upon the model predic¬ 
tion for the CGFF. 
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6.2.1 Galaxy mergers and interactions 

In GALFORM spheroids are created following galaxy mergers and 
disc instabilities. These events can also trigger starburst events, 
which would act to deplete the cold gas reservoirs of the merger 
remnant. We might expect therefore that varying the parameter gov¬ 
erning the timescale of these starbursts would have an effect upon 
the CGLF. We find, however, that changing the duration of star- 
bursts has little impact upon the predicted CGLF. 

From their analysis, Gonzalez et al. (2009) found that chang¬ 
ing the prescription for calculating the size of the stellar spheroid 
following a galaxy merger had a large impact upon the sizes of 
bright elliptical galaxies at low redshift. However, when we adopt 
their suggested parameter values we see little change in the CGLF 
or the predicted sizes for the most massive cluster galaxies. This 
might suggest that the amount of stellar mass produced in galaxy 
mergers is less important in high redshift cluster galaxies compared 
to the local Universe, which agrees with the lack of sensitivity we 
have seen to the starburst timescales. 

Font et al. (2008) demonstrated that the incorporation of grad¬ 
ual ram-pressure stripping into GALFORM improves the model pre¬ 
dictions for the colours of satellite galaxies compared with obser¬ 
vations. Recently, Lagos et al. (2014) have also shown that grad¬ 
ual ram pressure stripping is needed to reproduce the atomic and 
molecular gas contents of early-type galaxies. In GALFORM, when 
galaxies become satellites they have their hot gas stripped instan¬ 
taneously. Font et al. (2008) included a prescription to delay this 
stripping and allowed satellite galaxies to retain a fraction of their 
hot gas for a longer period, thus delaying the quenching of their 
star formation. However, when we include the Font et al. treatment 
for stripping, we again see a negligible change in the prediction for 
the CGLF, which might suggest that the galaxy sizes are having the 
greatest impact upon the dust attenuation calculation. 

In our reference model, the merger timescale for a satellite 
galaxy is calculated, based upon dynamical friction arguments, ev¬ 
ery time its host halo undergoes a merger event. The satellite is 
assumed to merge onto the central galaxy after this time, irrespec¬ 
tive of whether the sub-halo hosting the galaxy is still identifiable in 
the simulation. Campbell et al. (2015) recently showed that an al¬ 
ternative scheme, where the merger timescale is computed once the 
host sub-halo can no longer be identified, leads to a change in the 
model predictions for the stellar mass function at z = 0. We find 
that adopting this scheme does not improve the model predictions 
for the CGLF. 


6.2.2 Supernovae feedback 

The amount of cold gas in massive cluster galaxies could be low¬ 
ered by reducing the strength of feedback due to supernovae such 
that the galaxies undergo more star formation at earlier epochs, 
prior to them falling into the clusters. However, feedback due to 
supernovae is thought to affect the faint-end slope of the global 
galaxy luminosity function (e.g. Benson et al. 2003). As such, it is 
possible that reducing the strength of the supernova feedback will 
reduce the gas content of the galaxies, but will produce an undesir¬ 
able boost in the faint-end of the CGLF above the observations. 

In the upper two panels of Fig. 13 we show the effect of in¬ 
dependently varying the free parameters ahot and v^ot- Variation 
of either parameter clearly affects the normalisation of the CGLF. 
Understandably, the change in normalisation is greater at the faint- 
end. For the bright-end of the CGLF, there is little change in the 
normalisation as the parameters are decreased below their fiducial 
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Figure 13. Impact on H-band cluster galaxy luminosity function at z ^ 1.4 
as predicted by the GALFORM model when varying the parameters gov¬ 
erning SN and AGN feedback. The top panel shows the predictions when 
varying the SN feedback parameter «hot> the middle panel shows the pre¬ 
dictions when varying the SN feedback parameter Uhot an d the bottom 
panel shows the predictions when varying the AGN feedback parameter 
o C ooi- The parameters used by our reference model are indicated in the 
legend of each panel. All predictions assume the dust attenuation calcula¬ 
tion described in § 2.3. The halo mass and aperture used to select the semi- 
analytical cluster galaxies are indicated in each panel. The data points are 
the same as from the right-hand panel of Fig. 3. 
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values, suggesting that further change in the supernova feedback 
alone would have minimal impact. The values for the parameters 
could be reduced further beyond the range considered in Fig. 13, 
but this would cause an excess in the faint-end above the observa¬ 
tions. As such, we conclude that changing the strength of the su¬ 
pernova feedback alone is unable to fix the deficit around the break 
in the CGLF. 


6.2.3 AGN feedback 

Feedback due to active galactic nuclei (AGN) is expected to have 
a dramatic impact on galaxies residing in relatively massive sys¬ 
tems, like galaxy clusters, through the quenching of any active star 
formation. 

Bower et al. (2006) and Croton et al. (2006) were amongst 
the first to demonstrate that introducing feedback due to AGN into 
semi-analytic models helps reduce the number of bright galaxies in 
the models, thus improving the match to the bright-end of the ob¬ 
served global galaxy luminosity function. Additionally the action 
of AGN feedback leads to the models predicting a bi-modal colour- 
magnitude relation, similar to that observed in the SDSS (Gonzalez 
et al. 2009). Since the bright end of the global luminosity function 
is dominated by cluster galaxies, we expect that adjusting the pa¬ 
rameters controlling AGN feedback will affect the predicted CGLF. 

As discussed in §2.2, in GALFORM the strength of AGN feed¬ 
back is governed by the parameter a coo \. The effect of varying the 
value of a CO o\ on the cluster galaxy luminosity function is shown 
in the lower panel of Fig. 13 where we plot the H-band CGLF for 
a range of values of a coo u either side of the value of a C ooi = 0.6 
used in the reference model. As expected, we see that variations 
in the value of a coo i has a dramatic effect upon the sharpness of 
the break and normalisation of the bright-end of the CGLF. In con¬ 
trast, the change in the normalisation of the faint-end of the CGLF 
is negligible, especially when a coo i is reduced below the reference 
value. 

We can see from Fig. 13 that reducing a coo i brings the model 
predictions into better agreement with the observations, though 
even for a coo i <0.4 there is still a slight discrepancy around the 
break in the CGLF, suggesting that a further reduction in a coo i 
would be necessary. Examination of the CMR shows that adopt¬ 
ing a C ooi <0.4 has negligible impact upon the galaxy colours, 
with the slope and the zero-point of the RS (fitted to those galaxies 
with log 10 (SFR//i -1 M©yr -1 ) ^ —2) changing on the order of 
one per cent. However, as we shall see in §6.2.4, assuming such 
weak AGN feedback has a significant impact upon the model pre¬ 
dictions at z = 0, which indicates that solely reducing AGN feed¬ 
back is not an adequate solution and that some other mechanism 
must be changed, or introduced, if the model is to correctly predict 
the colours and abundances of high redshift cluster galaxies. 


6.2.4 Influence on local Universe predictions 

Following our brief parameter search we have found that the 
strength of AGN feedback could be used to reduce the discrepancy 
between the model predictions and observations of high redshift 
galaxy clusters. We now examine how our attempt to match ob¬ 
servations of high redshift clusters changes the predictions of the 
model at z = 0, which were originally used to calibrate the model. 

One of the principal statistics used to constrain the parameters 
of the GALFORM model is the global galaxy luminosity function 
(of both field and cluster galaxies) at z — 0, specifically in the 


bj and K-bands. In Fig. 14 we show the global galaxy luminosity 
function at z — 0 for the bj, z, H and K bands. The predictions for 
our reference GALFORM model (adopting the default dust attenua¬ 
tion calculation) are shown by the solid line. Since the parameters 
of the fiducial model have been constrained using the bj and K 
band luminosity functions, this model provides a good match to the 
global luminosity function in each of the four bands. 

The red, dashed line in Fig. 14 shows the impact at z = 0 
of assuming weaker dust attenuation, in this case a Calzetti et al. 
(2000) law with E(B — V) = 0.02, which is gives a better match 
to the CGLF. As expected, weaker dust attenuation boosts the abun¬ 
dance of bright galaxies, with a reduction in the predicted counts 
around the knee of the luminosity function. The impact on the lumi¬ 
nosity function of adopting weaker dust attenuation becomes more 
significant as one moves from the near-infrared towards the opti¬ 
cal, with the K-band showing the smallest change and the bj-band 
showing the largest. We have seen in §6.2.3 that weaker AGN feed¬ 
back is necessary to reconcile the model predictions for the CGLF 
and the observations. As such, we also plot in Fig. 14 the predicted 
z — 0 luminosity function for the model when a coo i = 0.3 is 
assumed. The effect on the luminosity function is dramatic, with 
the weaker AGN feedback leading to a significant excess of bright 
galaxies. 

These latter two predictions again highlight the challenge fac¬ 
ing current galaxy formation models. In addition to reproducing the 
sizes, luminosities and colours of massive galaxies in clusters, the 
models need to be able to remain consistent with observations of 
the local Universe. To achieve this requires more than one parame¬ 
ter to be varied as well as the possible inclusion of new physics. 


7 SUMMARY & CONCLUSIONS 

We have compiled observations of high redshift (z > 1) galaxy 
clusters, which we compare to the predictions of the Gonzalez- 
Perez et al. (2014) variant of the GALFORM semi-analytical galaxy 
formation model, which we treat as our reference model. The 
statistics that we consider are the cluster galaxy luminosity func¬ 
tion (CGLF) and the colour-magnitude relation (CMR). To identify 
cluster galaxies in the semi-analytic catalogue, we select only those 
galaxies in halos with mass greater than 1.2 x 10 14 /i _ 1 M©. We 
further use the distant observer approximation to apply an aperture 
and reject those galaxies lying further than 120" away from the 
halo centre. 

Our reference GALFORM model predicts a CGLF in reason¬ 
able agreement with the observed CGLF at the faint and bright 
ends, but significantly under-predicts the number of cluster galax¬ 
ies around the break in the CGLF. Examination of several possible 
factors that might affect the model predictions, including aperture 
size and halo mass selection, indicates that the discrepancy between 
the observations and the model predictions is likely caused by the 
reference model applying an overly large dust attenuation. If we 
instead apply a weaker dust attenuation, which we represent using 
a Calzetti et al. (2000) law, then the reference model prediction is 
able to provide a much better fit to the observed CGLF. We note that 
we do not advocate that a Calzetti et al. law is the correct descrip¬ 
tion for dust attenuation at high redshift, but instead have simply 
used the law to demonstrate the impact of the large dust attenuation 
predicted in our reference model. 

In contrast, the reference GALFORM model predicts a CMR 
that is qualitatively consistent with the observed colours of clus¬ 
ter galaxies at z ~ 1.2, z ~ 1.4 and z ~ 1.6. We provide lin- 
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Figure 14. Predicted luminosity functions at z = 0 for the global galaxy population (both field and cluster galaxies). The panels show the luminosity functions 
in the bj, z, H and K-bands, as labelled. Data points show the observational estimates from the Two-Micron All-Sky Survey (2MASS, Kochanek et al. 2001), 
Two-degree Field Galaxy Redshift Survey 2dFGRS (2dFGRS, Cole et al. 2001; Norberg et al. 2002), the Six-degree Field Galaxy Survey (6dFGS, Jones et al. 
2006) and the Galaxy And Mass Assembly Survey (GAMA, Driver et al. 2012). The solid blue line shows the prediction for the fiducial GALFORM model 
(Gonzalez-Perez et al. 2014) with the default dust attenuation calculations. The dashed line shows the prediction of the fiducial model when, instead of the 
default extinction, a Calzetti et al. (2000) law with E(B — V) = 0.02 is adopted. The green, dotted line shows the GALFORM prediction, with the default dust 
extinction, when the parameter a coo i is reduced from 0.6 to 0.3. 


ear fits to the red sequence, RS, for different colour-spaces using a 
subset of passively evolving galaxies selected using the sSFR cut, 
log 10 (sSFR/Gyr -1 ) ^ —1. The slopes and zero-points of these 
fits are broadly consistent with observationally derived estimates. 

The CMR predicted by our reference model displays a subset 
of very red galaxies, which appear in a ‘plume’ above the predicted 
RS. We determine that these galaxies are a result of the large dust 
attenuation that is causing the discrepancy between the observed 
CGLF and that in the model. However, although a weaker dust at¬ 


tenuation improves the model predictions for the CGLF, assuming 
a weaker dust attenuation worsens the predicted RS, with the cre¬ 
ation of a branch of blue galaxies extending below the RS. 

Examination of the properties of those galaxies with large 
dust attenuation reveals them to be highly star-forming, with large 
amounts of cold gas and with small scale sizes. The large reser¬ 
voirs of cold gas and the compact sizes appear to be the cause of 
the large predicted dust attenuation. To gain some insight into this 
problem we briefly examined how varying several key model pa- 
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rameters changes the predicted cluster statistics, in particular the 
CGLF. We find that a reduction in the strength of feedback due 
to AGN is able to provide some improvement in the form of the 
CGLF, but makes the colours of the model galaxies too blue. In 
addition, introducing a weaker AGN feedback significantly affects 
the model predictions at z = 0 by boosting the number of bright 
galaxies and over-predicting the counts at the bright end of z — 0 
global galaxy luminosity function. 

There are several possible explanations for the discrepancy 
between the observations and the model predictions. Amongst the 
most probable are that GALFORM is under-predicting the star for¬ 
mation of these massive cluster galaxies, or allowing too much gas 
cooling, which would leave them too faint (and most likely too 
blue) and with lots of cold gas by the time they become cluster 
galaxies. The problem of semi-analytical models under-predicting 
star formation at high redshift has been commented on several 
times in the literature (e.g. Daddi et al. 2007; Damen et al. 2009; 
Fontanot et al. 2009a; Dutton et al. 2010; Weinmann et al. 2011, 
2012). Similar deficiencies have also been reported in hydrody- 
namical simulations (e.g. Kannan et al. 2014; Wang et al. 2015). 
However, it is also worth noting a recent result from Chang et al. 
(2015) who suggest that observational calibration of SFR estimates 
could be wrong by approximately a factor of two, which would re¬ 
sult in previous observational SFR estimates being a factor of two 
too large. If this is the case then correcting for this would bring the 
observational SFR estimates and theoretical predictions into closer 
agreement. Additionally, the under-prediction of the sizes of the 
galaxies in GALFORM is likely to be having a significant impact 
upon the other model predictions. Predicting correct galaxy sizes is 
a long standing problem for semi-analytical models since it is en¬ 
tangled with the contraction of the host halo. These conclusions are 
in agreement with previous comparisons to semi-analytical model 
predictions (e.g. Gonzalez et al. 2009; Weinzirl et al. 2014). 

Overall, these results demonstrate the challenge facing cur¬ 
rent galaxy formation models such as GALFORM: how to match the 
luminosity abundance and colours of cluster galaxies, whilst re¬ 
maining consistent with the observed properties of galaxies in the 
local Universe. Achieving this will require incorporating into the 
models a better understanding of galaxy evolution in extreme envi¬ 
ronments. However, we must keep in mind that the comparisons in 
this work are limited to a small sample of individual galaxy clus¬ 
ters and that understanding of the astrophysical processes affecting 
the high redshift galaxy population remains uncertain. Ultimately, 
our understanding of the abundance and properties of high redshift 
clusters will only improve as we improve our statistics with up and 
coming deep, wide-field galaxy surveys such as DES and Euclid. 
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